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The running coupling constant can be estimated by computing gluon two- and three-point 
Green functions from the lattice. Computing in lattice implies working in a fixed gauge sector 
(Landau). An source of systematic uncertainty is then the contribution of non gauge- invariant 
condensates, as (A 2 ), generating power corrections that might still be important at very high 
energies. We study the impact of this gluon condensate on the analysis of gluon propagator 
and vertex lattice data and on the estimate of Aj-jg. We finally try a qualitative description 
of this gluon condensate through the instanton picture. 



1 The running coupling constant from the lattice 

That the running coupling constant can be extracted from the three-gluon vertex in the Landau 
gauge was proposed several years ago in a seminal worklll. The key lies on the appropriate choice 
of renormalisation scheme: the so-called Momentum Substraction (MOM) schemes are defined 
such that all the renormalised Green functions take their tree-level expressions after replacing 
bare by renormalised constants. Then 



9R{k ] ~ (Gvmf ( z {k ] i (1) 

where G^ and G® are the scalar form factor of two- and three-gluon Green functions and Z 
is the gluon propagator renormalisation constant. Alternative kinematics for the renormalisation 
point are possible, mainly two among them: symmetric (k 2 = k\ = k\ = k 2 ) and asymmetric 
(k 2 = k 2 = k 2 ,k 2 = 0). Analysis of lattice computations of these Green functions renormalised 



in both yielded □ the estimates of the coupling and A^ collected in Table ||. Concerned as 
we were by the determination and renormalisation of the gluon propagator for the analysis of 
the three-point Green Function, we exploited @ the first by matching the data to the three- 
loop perturbative prediction and hence estimate a purely perturbative coupling and A^g. The 
disagreement between estimates from both gluon vertex and propagator (see Table [l]) manifests 
the impact of some ucontrolled systematic uncertainty. Even worse, the so estimated A^g does 
not behave as a scale invariant!! 

Table 1: Estimates of the coupling and A^r from two- and three-point Green function methods . 



a(4.3GeV) 


Three-point 


Two-point 


0.269(3) 




A^(4.3GeV) 


299(7) 




a(9.6GeV) 


0.176(2) 


0.193(3) 


A^(9.6GeV) 


266(7) 


319(14) 



If one empirically adds a 1/p 2 power correction to the perturbative formulae (we work at 
three loops) 

G(p 2 ) = G(p 2 ) 3 ioop S (l + ^) > a ^P 2 ) = "Whoops (l + , (2) 

success is gained both in obtaining estimates by two- and three-point methods that agree to 
each other and in restoring the scale invariance of the estimated Aj^g parameter! (see fig. |j].a). 
A third outcome arises from matching to Eq. [||: the estimate of A-^g results 237(4) MeV 
(only the statistical error is quoted), in astonishing agreement with Schrdedinger functional' 'si: 
238(19) MeV !! 




Figure 1: In plot (a), the upper solid points are for the estimates of obtained through the non perturbative 
couplings directly computed from lattice Green functions. The lower ones are obtained through aj crt in Eq. |^ 
by the substraction of the fitted power contribution. Plot (b) shows the best fit of lattice data for the coupling 

renormalised in asymmetric MOM to Eqs. B, H. 
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Figure 2: Diagrams in the OPE expansion of gluon two- and three-point Green functions. 
2 The O.P.E. picture 

How large should c and d are? The Operator Product Expansion (OPE) can be invoked to let us 
gain some physical insight into the matter of these coefficients. The sum rule approach connects 
the power corrections to any non-local matrix element and the QCD vacuum expectation value of 
the local operators from the expansion, the so-called condensates, and gives the prescription to 
compute perturbatively the coefficients of the expansion (Wilson coefficients). Symmetries and 
power counting tell us that the gluon condensate (A 2 ) gives the major contribution to our gluon 
Green functions. This is of course a non gauge-invariant contribution but, as far as we work in 
a fixed gauge, Landau gauge, non-vanishing contributions from gauge-dependant condensates, 
i.e. (A 2 ), should be expectedi. Then, our coefficients c and d can be computedB (see fig. ||| ) 
and written in terms of the gluon condensate renormalised at any momentum scale [/,, 



TO +70 ^ , _ 7Q + 7Q 



with 70 = 13/2, fa = ll,7o = 3JVc/4. The factor x depends on the particular kinematics we 
choose for the vertex renormalisation point; case asymmetric: x = 1, case symmetric: x = 3. 

Then, Eqs. [2L |3| are at our disposal for trying a coherent description of gluon propagator 
and vertex lattice data, with the sole necessity of fitting two ingredients of non-perturbative 
nature: (A 2 ) and A-^g (except for an overall factor in the propagator analysis). We expand 
purely perturbative series up to the third loop and only keep the leading logarithm contribution 
for the Wilson coefficient. Of course this does not guarantee us to reach the asymptotic regime 
where expansions well behave. Thus, we approach the problem through a combined analysis of 
both gluon propagator and vertex, where we look for fitting the same Aj^g parameter in both 
and two independent estimates of {A 2 ) to be compared. The results of applying this to our 
lattice dataQ are in Table 0. The best fit parameters lead, for instance in the case of the MOM 
asymmetric vertex, to the plot (b) in fig. EL 

Consequently, the contribution of this non gauge-invariant condensate (A 2 ) seems to explain 
rather well the systematic deviation in the matchings of our gluon propagator and vertex lattice 
data to purely perturbative formulae. 



Table 2: Comparison of results obtained to the renormalisation momentum scale fi = 10 GeV. 





asymmetric MOM 


symmetric MOM 




260(18) MeV 


233(28) MeV 


WW)**} 

L V J prop 


1.39(14) GeV 


1.55(17) GeV 




2.3(6) GeV 


1.9(3)GeV 



3 Instantons and (A 2 ) condensate 

The physical origin of this gluon condensate is a major question. A common belief is that 
an instanton ensemble (liquid or gas) provides with a fair description of important features of 
the QCD vacuum. Then, if one considers a hard gluon of momentum p^ propagating in an 
instanton gas background (we crudely assume that instantons do not interact to each other), 
the propagator can be computed with Feynman graphs and it is easy to see that the dominant 
contribution coming from its interaction wiht the instanton gauge field is 0(l/p 2 ) when the gauge 
field momentum, kn « p^- This correction is equal to the standard OPE Wilson coefficient for 
the propagator times 

CO = / = i^V^± (4) 

where A" is the standard 't Hooft-Polyakov solution in the singular Landau eauge, nj(nA) is 
the number of (anti-)instantons and p is the instanton radius. We estimates this instanton- 
induced condensate to be 1.76(23) GeV 2 by performing several simulations at (3 = 6.0 on a 24 4 
lattice, applying the cooling procedure to them for killing UV fluctuations and computing the 
instanton density to be applied in Eq. ||. This result is to be compared with the OPE estimate 
from the propagator analysis. We should first run down ]] the value in Table Q to the lowest 
renormalisation momentum scale we still believe for our OPE analysis, that is ~ 2.5 GeV. We 
obtain 1.4(3) GeV 2 . This fair agreement seems to indicate that the gluon condensate receives a 
significant instantonic contribution. 
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a In the instanton analysis, the renormalisation is performed through the cooling that kills the UV fluctuations. 
The only new physical momentum scale related to the instanton background emerging after the cooling procedure 
is the inverse of the instanton radius ~ 0.5 GeV 



